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Cross-species transmission of viruses from wildlife animal 
reservoirs poses a marked threat to human and animal health!. 
Bats have been recognized as one of the most important reservoirs 
for emerging viruses and the transmission of a coronavirus 
that originated in bats to humans via intermediate hosts was 
responsible for the high-impact emerging zoonosis, severe acute 
respiratory syndrome (SARS)^!?, Here we provide virological, 
epidemiological, evolutionary and experimental evidence that 
a novel HKU2-related bat coronavirus, swine acute diarrhoea 
syndrome coronavirus (SADS-CoV ), is the aetiological agent that 
was responsible for a large-scale outbreak of fatal disease in pigs in 
China that has caused the death of 24,693 piglets across four farms. 
Notably, the outbreak began in Guangdong province in the vicinity 
of the origin of the SARS pandemic. Furthermore, we identified 
SADS-related CoVs with 96-9896 sequence identity in 9.896 (58 out 
of 591) of anal swabs collected from bats in Guangdong province 
during 2013-2016, predominantly in horseshoe bats (Rhinolophus 
spp.) that are known reservoirs of SARS-related CoVs. We found 
that there were striking similarities between the SADS and SARS 
outbreaks in geographical, temporal, ecological and aetiological 
settings. This study highlights the importance of identifying 
coronavirus diversity and distribution in bats to mitigate future 
outbreaks that could threaten livestock, public health and economic 
growth. 

The emergence of SARS in southern China in 2002, which was 
caused by a previously unknown coronavirus (SARS-CoV)! and 
has led to more than 8,000 human infections and 774 deaths (http:// 
www.who.int/csr/sars/en/), highlights two new frontiers in emerging 
infectious diseases. First, it demonstrates that coronaviruses are capable 
of causing fatal diseases in humans. Second, the identification of bats as 
the reservoir for SARS-related coronaviruses, and the fact that SARS- 
CoV?-? probably originated in bats, firmly establishes that bats are an 
important source of highly lethal zoonotic viruses, such as Hendra, 
Nipah, Ebola and Marburg viruses!*. 

Here we report on a series of fatal swine disease outbreaks in 
Guangdong province, China, approximately 100 km from the location 
of the purported index case of SARS. Most strikingly, we found that 
the causative agent of this swine acute diarrhoea syndrome (SADS) is 
a novel HKU2-related coronavirus that is 98.48% identical in genome 
sequence to a bat coronavirus, which we detected in 2016 in bats in a 
cave in the vicinity of the index pig farm. This new virus (SADS-CoV) 


originated from the same genus of horseshoe bats (Rhinolophus) as 
SARS-CoV. 

From 28 October 2016 onwards, a fatal swine disease outbreak was 
observed in a pig farm in Qingyuan, Guangdong province, China, 
very close to the location of the first known index case of SARS in 
2002, who lived in Foshan (Extended Data Fig. 1a). Porcine epidemic 
diarrhoea virus (PEDV, a coronavirus) had caused prior outbreaks at 
this farm, and was detected in the intestines of deceased piglets at the 
start of the outbreak. However, PEDV could no longer be detected in 
deceased piglets after 12 January 2017, despite accelerating mortality 
(Fig. 1a), and extensive testing for other common swine viruses yielded 
no results (Extended Data Table 1). These findings suggested that this 
was an outbreak of a novel disease. Clinical signs are similar to those 
caused by other known swine enteric coronaviruses!” !? and include 
severe and acute diarrhoea and acute vomiting, leading to death due 
to rapid weight loss in newborn piglets that are less than five days of 
age. Infected piglets died 2-6 days after disease onset, whereas infected 
sows suffered only mild diarrhoea and most sows recovered within 
two days. The disease caused no signs of febrile illness in piglets or 
sows. The mortality rate was as high as 9096 in piglets that were five 
days or younger, whereas in piglets that were older than eight days, the 
mortality dropped to 596. Subsequently, SADS-related outbreaks were 
found in three additional pig farms within 20-150 km ofthe index farm 
(Extended Data Fig. 1a) and, by 2 May 2017, the disease had caused 
the death of 24,693 piglets at these four farms (Fig. 1a). In farm A 
alone, 6496 (4,659 out of 7,268) of all piglets that were born in February 
died. The outbreak has abated, and measures that were taken to control 
SADS included separation of sick sows and piglets from the rest of the 
herd. A qPCR test described below was used as the main diagnostic 
tool to confirm SADS-CoV infection. 

A sample collected from the small intestine of a diseased piglet was 
analysed by metagenomics analysis using next-generation sequenc- 
ing (NGS) to identify potential aetiological agents. Of the 15,256,565 
total reads obtained, 4,225 matched sequences of the bat CoV HKU2, 
which was first detected in Chinese horseshoe bats in Hong Kong and 
Guangdong province, China. By de novo assembly and targeted PCR, 
we obtained a 27,173-bp CoV genome that shared 9596 sequence iden- 
tity to HKU2-CoV (GenBank accession number NC. 009988). Thirty- 
three full genome sequences of SADS-CoV were subsequently obtained 
(8 from farm A, 5 from farm B, 11 from farm C and 9 from farm D) that 
were 99.996 identical to each other (Supplementary Table 1). 
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Fig. 1 | Detection of SADS-CoV infection in pigs in Guangdong, China. 
a, Records of daily death toll on the four farms from 28 October 2016 to 

2 May 2017. b, Detection of SADS-CoV by qPCR. The y axis shows the 
log(copy number per 10° copies of 18S rRNA). n = 12 sick piglets, 5 sick 
sows, 16 recovered sows and 10 healthy piglets. c, Tissue distribution of 
SADS-CoV in diseased pigs. n = 3. Data are mean c s.d.; dots represent 


Using qPCR targeting the nucleocapsid gene (Supplementary 
Table 2), we detected SADS-CoV in acutely sick piglets and sows, but 
not in recovered or healthy pigs on the four farms, nor in nearby farms 
that showed no evidence of SADS. The virus replicated to higher titres 
in piglets than in sows (Fig. 1b). SADS-CoV displayed tissue tropism 
of the small intestine (Fig. 1c), as observed for other swine enteric 
coronaviruses”. Retrospective PCR analysis revealed that SADS-CoV 
was present on farm A during the PEDV epidemic, where the first 
strongly positive SADS-CoV sample was detected on 6 December 2016. 
From mid-January onwards, SADS-CoV was the dominant viral agent 
detected in diseased animals (Extended Data Fig. 1b). It is possible 
that the presence of PEDV early in the SADS-CoV outbreak may have 
somehow facilitated or enhanced spillover and amplification of SADS. 
However the fact that the vast majority of piglet mortality occurred 
after PEDV infection had become undetectable suggests that SADS- 
CoV itself causes a lethal infection in pigs that was responsible for these 
large-scale outbreaks, and that PEDV does not directly contribute to 
its severity in individual pigs. This was supported by the absence of 
PEDV and other known swine diarrhoea viruses during the peak and 
later phases of the SADS outbreaks in the four farms (Extended Data 
Table 1). 

We rapidly developed an antibody assay based on the S1 domain of 
the spike (S) protein using a luciferase immunoprecipitation system?!. 
Because SADS occurs acutely and has a rapid onset in piglets, serolog- 
ical investigation was conducted only in sows. Among 46 recovered 
sows tested, 12 were seropositive for SADS-CoV within three weeks 
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individual values. d, Detection of SADS-CoV antibodies. n = 46 sows 
from whom serum was first taken in the first three weeks of the outbreak 
(First bleed), n = 8 sows from whom serum was taken again (Second 
bleed) at more than one month after the onset of the outbreak, n = 8 sera 
from healthy pig controls, n — 35 human sera from pig farmers. 


of infection (Fig. 1d). To investigate possible zoonotic transmission, 
serum samples from 35 farm workers who had close contact with sick 
pigs were also analysed using the same luciferase immunoprecipitation 
system approach and none were positive for SADS-CoV. 

Although the overall genome identity of SADS-CoV and HKU2- 
CoV is 95%, the S gene sequence identity is only 86%, suggesting that 
the previously reported HKU2-CoV is not the direct progenitor of 
SADS-CoV, but that they may have originated from a common ances- 
tor. To test this hypothesis, we developed a SADS-CoV-specific qPCR 
assay based on its RNA-dependent RNA polymerase (RdRp) gene 
(Supplementary Table 2) and screened 591 bat anal swabs collected 
between 2013 and 2016 from seven different locations in Guangdong 
province (Extended Data Fig. 1a). A total of 58 samples (9.896) tested 
positive (Extended Data Table 2), all were from Rhinolophus spp. bats 
that are also the natural reservoir hosts of SARS-related coronavi- 
ruses’ ?, Four complete genome sequences with the highest RdRp 
PCR-fragment sequence identity to that of SADS-CoV were deter- 
mined by NGS. They are very similar in size (27.2 kb) compared to 
SADS-CoV (Fig. 2a) and we tentatively call them SADS-related corona- 
viruses (SADSr-CoV ). Overall sequence identity between SADSr-CoV 
and SADS-CoV ranges from 96 to 98%. Most importantly, the S protein 
of SADS-CoV shared more than 98% sequence identity with sequences 
of two of the SADSr-CoVs (samples 162149 and 141388), compared to 
86% with HKU2-CoV. The major sequence differences among the four 
SADSr-CoV genomes were found in the predicted coding regions of 
the S and NS7a and NS7b genes (Fig. 2a). In addition, the coding region 
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162140-CoV: 27,177 bp, 98.4896 identity 


141388-CoV: 27,174 bp, 98.0596 identity 


8462-CoV: 27,200 bp, 96.36% identity 


8495-CoV: 27,198 bp, 96.28% identity 


HKU2-CoV-GD: 27,165 bp, 95.0996 identity 


Fig. 2 | Genome and phylogenetic analysis of SADS-CoV and SADSr- 
CoV. a, Genome organization and comparison. Colour-coding for different 
genomic regions as follows. Green, non-structural polyproteins ORF1a and 
ORFIb; yellow, structural proteins S, E, M and N; blue, accessory proteins 
NS3a, NS7a and NS7b; Orange, untranslated regions. The level of sequence 
identity of SADSr-CoV to SADS-CoV is illustrated by different patterns 


of the S protein N-terminal (S1) domain was determined from 19 bat 
SADSr-CoVs to enable more detailed phylogenetic analysis. 

The phylogeny of S1 and the full-length genome revealed a high 
genetic diversity of alphacoronaviruses among bats and strong coev- 
olutionary relationships with their hosts (Fig. 2b and Extended Data 
Fig. 2), and showed that SADS-CoVs were more closely related to 
SADSr-CoVs from Rhinolophus affinis than from Rhinolophus sinicus, 
in which HKU2-CoV was found. Both phylogenetic and haplotype net- 
work analyses demonstrated that the viruses from the four farms proba- 
bly originated from their reservoir hosts independently (Extended Data 
Fig. 3), and that a few viruses might have undergone further genetic 
recombination (Extended Data Fig. 4). However, molecular clock 
analysis of the 33 SADS-CoV genome sequences failed to establish a 
positive association between sequence divergence and sampling date. 
Therefore, we speculate that either the virus was introduced into pigs 
from bats multiple times, or that the virus was introduced into pigs 
once, but subsequent genetic recombination disturbed the molecular 
clock. 

For viral isolation, we tried to culture the virus in a variety of cell lines 
(see Methods for details) using intestinal tissue homogenates as starting 
material. Cytopathogenic effects were observed in Vero cells only after 
five passages (Extended Data Fig. 5a, b). The identity of SADS-CoV was 
verified in Vero cells by immunofluorescence microscopy (Extended 
Data Fig. 5c, d) and by whole-genome sequencing (GenBank accession 
number MG557844). Similar results were obtained by other groups’? ^. 

Known coronavirus host cell receptors include angiotensin-con- 
verting enzyme 2 (ACE2) for SARS-related CoV, aminopeptidase N 
(APN) for certain alphacoronaviruses, such as human (H)CoV-229E, 
and dipeptidyl peptidase 4 (DPP4) for Middle East respiratory syn- 
drome (MERS)-CoV**™®. To investigate the receptor usage of SADS- 
CoV, we tested live or pseudotyped SADS-CoV infection on HeLa 
cells that expressed each of the three molecules. Whereas the positive 
control worked for SARS-related CoV and MERS-CoV pseudoviruses, 
we found no evidence of enhanced infection or entry for SADS-CoV, 
suggesting that none of these receptors functions as a receptor for virus 
entry for SADS-CoV (Extended Data Table 3). 

To fulfill Kochs postulates for SADS-CoV, two different types of 
animal challenge experiments were conducted (see Methods for 
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details). The first challenge experiment was conducted with specific 
pathogen-free piglets that were infected with a tissue homogenate of 
SADS-CoV-positive intestines. Two days after infection, 3 out of 7 
animals died in the challenge group whereas 4 out of 5 survived in 
the control group. Incidentally, the one piglet that died in the control 
group was the only individual that did not receive colostrum due to a 
shortage in the supply. It is thus highly likely that lack of nursing and 





Fig. 3 | Inmunohistopathology of SADS-CoV infected tissues. a- 

d, Sections of jejunum tissue from control (a, c) and infected (b, d) farm 
piglets four days after inoculation were stained with haematoxylin and 
eosin (a, b) or rabbit anti-SADSr-CoV N serum (red), DAPI (blue) and 
mouse antibodies against epithelial cell markers cytokeratin 8, 18 and 19 
(green) in (c, d). SADS-CoV N protein is evident in epithelial cells and 
deeper in the tissue of infected piglets, which exhibit villus shortening. 
Scale bars, 200 um (a, b) and 50 jum (c, d). The experiment was conducted 
three times independently with similar results. 
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inability to access colostrum was responsible for the death (Extended 
Data Table 4). For the second challenge, healthy piglets were acquired 
from a farm in Guangdong that had been free of diarrheal disease for 
a number of weeks before the experiment, and were infected with the 
cultured isolate of SADS-CoV or tissue-culture medium as control. 
Of those inoculated with SADS-CoV, 5096 (3 out of 6) died between 
2 and 4 days after infection, whereas all control animals survived 
(Extended Data Table 5). All animals in the infected group suffered 
watery diarrhoea, rapid weight loss and intestinal lesions (determined 
after euthanasia upon experiment termination, Extended Data Tables 4, 
5). Histopathological examination revealed marked villus atrophy in 
SADS-CoV inoculated farm piglets four days after inoculation but not 
in control piglets (Fig. 3a, b) and viral N protein-specific staining was 
observed mainly in small intestine epithelial cells of the inoculated 
piglets (Fig. 3c, d). 

The current study highlights the value of proactive viral discovery 
in wildlife, and targeted surveillance in response to an emerging infec- 
tious disease event, as well as the disproportionate importance of bats 
as reservoirs of viruses that threaten veterinary and public health’. It 
also demonstrates that by using modern technological platforms, such 
as NGS, luciferase immunoprecipitation system serology and phyloge- 
netic analysis, key experiments that traditionally rely on the isolation of 
live virus can be performed rapidly before virus isolation. 








Online content 

Any Methods, including any statements of data availability and Nature Research 
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METHODS 


Sample collection. Bats were captured and sampled in their natural habitat in 
Guangdong province (Extended Data Fig. 1) as described previously*. Faecal swab 
samples were collected in viral transport medium (V TM) composed of Hanks 
balanced salt solution at pH 7.4 containing BSA (1%), amphotericin (15 ug mI P), 
penicillin G (100 units ml~') and streptomycin (50 ug ml). Stool samples from 
sick pigs were collected in VTM. When appropriate and feasible, intestinal samples 
were also taken from deceased animals. Samples were aliquoted and stored at -80 
°C until use. Blood samples were collected from recovered sows and workers on the 
farms who had close contact with sick pigs. Serum was separated by centrifugation 
at 3,000g for 15 min within 24 h of collection and preserved at 4 °C. Human serum 
collection was approved by the Medical Ethics Committee ofthe Wuhan School of 
Public Health, Wuhan University and Hummingbird IRB. Human, pigs and bats 
were sampled without gender or age preference unless indicated (for example, 
piglets or sows). No statistical methods were used to predetermine sample size. 
Virus isolation. The following cells were used for virus isolation in this study: Vero 
(cultured in DMEM and 1096 FBS); Rhinolophus sinicus primary or immortalized 
cells generated in our laboratory (all cultured in DMEM/F12 and 15% FBS): kidney 
primary cells (RsKi9409), lung primary cells (RsLu4323), lung immortalized cells 
(RsLuT), brain immortalized cells (RsBrT) and heart immortalized cells (RsHeT); 
and swine cell lines: two intestinal porcine enterocytes cell lines, IPEC (RPMI1640 
and 1096 FBS) and SIEC (DMEM and 1096 FBS), three kidney cell lines PK15, 
LLC-PK1 (DMEM and 1096 FBS for both) and IBRS (MEM and 1096 FBS), and 
one pig testes cell line, ST (DMEM and 10% FBS). All cell lines were tested free of 
mycoplasma contamination, species were confirmed and authenticated by micro- 
scopic morphologic evaluation. None of the cell lines was on the list of commonly 
misidentified cell lines (by the ICLAC). 

Cultured cell monolayers were maintained in their respective medium. PCR- 

positive pig faecal samples or the supernatant from homogenized pig intestine 
(in 200 ul VTM) were spun at 8,000g for 15 min, filtered and diluted 1:2 with 
DMEM supplemented with 16 ug ml! trypsin before addition to the cells. After 
incubation at 37 °C for 1 h, the inoculum was removed and replaced with fresh 
culture medium containing antibiotics (below) and 16 pg ml“! trypsin. The cells 
were incubated at 37 °C and observed daily for cytopathic effect (CPE). Four blind 
passages (three-day interval between every passage) were performed for each sam- 
ple. After each passage, both the culture supernatant and cell pellet were exam- 
ined for the presence of virus by RT-PCR using the SADS-CoV primers listed in 
Supplementary Table 2. Penicillin (100 units ml~*) and streptomycin (15 ug ml~’) 
were included in all tissue culture media. 
RNA extraction, S1 gene amplification and qPCR. Whenever commercial kits 
were used, the manufacturer's instructions were followed without modification. 
RNA was extracted from 200 ul of swab samples (bat), faeces or homogenized 
intestine (pig) with the High Pure Viral RNA Kit (Roche). RNA was eluted in 50 
ul of elution buffer and used as the template for RT-PCR. Reverse transcription 
was performed using the SuperScript III kit (Thermo Fisher Scientific). 

To amplify S1 genes from bat samples, nested PCR was performed with prim- 
ers designed based on HKU2-CoV (GenBank accession number NC, 009988.1)!? 
(Supplementary Table 2). The 25-l first-round PCR mixture contained 2.5 ul 
10x PCR reaction buffer, 5 pmol of each primer, 50 mM MgCh, 0.5 mM dNTP, 
0.1 ul Platinum Taq Enzyme (Thermo Fisher Scientific) and 1 ul cDNA. The 50-l 
second-round PCR mixture was identical to the first-round PCR mixture except 
for the primers. Amplification of both rounds was performed as follows: 94 °C for 
5 min followed by 60 cycles at 94 °C for 30 s, 50 °C for 40 s, 72 °C for 2.5 min, and a 
final extension at 72 °C for 10 min. PCR products were gel-purified and sequenced. 

For qPCR analysis, primers based on SADS-CoV RdRp and N genes were used 
(Supplementary Table 2). RNA extracted from above was reverse-transcribed using 
PrimeScript RT Master Mix (Takara). The 10 ul qPCR reaction mix contained 5 
ul 2x SYBR premix Ex Tagll (Takara), 0.4 uM of each primer and 1 ul cDNA. 
Amplification was performed as follows: 95 °C for 30 s followed by 40 cycles at 95 
°C for 5 s, 60 °C for 30 s, and a melting curve step. 

Luciferase immunoprecipitation system assay. The SADS-CoV S1 gene was 
codon-optimized for eukaryotic expression, synthesized (GenScript) and cloned 
in frame with the Renilla luciferase gene (Rluc) and a Flag tag in the pREN2 vec- 
tor’. pREN2-S1 plasmids were transfected into Cos-1 cells using Lipofectamine 
2000 (Thermo Fisher Scientific). At 48 h post-transfection, cells were collected, 
lysed and a luciferase assay was performed to determine Rluc expression for both 
the empty vector (pREN2) and the pREN2-S1 construct. For testing of unknown 
pig or human serum samples, 1 ul of serum was incubated with 10 million units of 
Rluc alone (vector) or Rluc-S1, respectively, together with 3.5 ul of a 3096 protein 
A/G UltraLink resin suspension (Pierce, Thermo Fisher Scientific). After extensive 
washing to remove unbounded luciferase-tagged antigens, the captured luciferase 
amount was determined using the commercial luciferase substrate kit (Promega). 
The ratio of Rluc-S1:Rluc (vector) was used to determine the specific S1 reactivity 
of pig and human sera. Commercial Flag antibody (Thermo Fisher Scientific) 
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was used as the positive control, and various pig sera (from uninfected animals 
in China or Singapore; or pigs infected with PEDV, TGEV or Nipah virus) were 
used as a negative control. 

Protein expression and antibody production. The N gene from SADSr-CoV 3755 
(GenBank accession number MF094702), which shares a 98% amino acid sequence 
identity to the SADS-CoV N protein, was inserted into pET-28a+ (Novagen) for 
prokaryotic expression. Transformed Escherichia coli were grown at 37 °C for 12-18 
h in medium containing 1 mM IPTG. Bacteria were collected by centrifugation 
and resuspended in 30 ml of 5 mM imidazole and lysed by sonication. The lysate, 
from which N protein expression was confirmed with an anti-His-tag antibody, 
was applied to Ni?” resin (Thermo Fisher Scientific). The purified N protein, at a 
concentration of 400 ig ml~', was used to immunize rabbits for antibody produc- 
tion following published methods”. After immunization and two boosts, rabbits 
were euthanized and sera were collected. Rabbit anti-N protein serum was used 
1:10,000 for subsequent western blots. 

Amplification, cloning and expression of human and swine genes. Construction 
of expression clones for human ACE2 in pcDNA3.1 has been described previously” 
28. Human DPP4 was amplified from human cell lines. Human APN (also known 
as ANPEP) was commercially synthesized. Swine APN (also known as ANPEP), 
DPP4 and ACE2 were amplified from piglet intestine. Full-length gene fragments 
were amplified using specific primers (provided upon request). Human ACE2 was 
cloned into pCDNA3.1 fused with a His tag. Human APN and DPP4, swine APN, 
DPP4 and ACE2 were cloned into pCAGGS fused with an S tag. Purified plasmids 
were transfected into HeLa cells. After 24 h, expression human or swine genes in 
HeLa cells was confirmed by immunofluorescence assay using mouse anti-His 
tag or mouse anti-S tag monoclonal antibodies (produced in house) followed by 
Cy3-labelled goat anti-mouse/rabbit IgG (Proteintech Group). 

Pseudovirus preparation. The codon-humanized S genes of SADS-CoV or MERS- 
CoV cloned into pcDNA3.1 were used for pseudovirus construction as described 
previously? ?*. In brief, 15 ug of each pHIV-Luc plasmid (pNL4.3.Luc.R-E-Luc) 
and the S-protein-expressing plasmid (or empty vector control) were co-trans- 
fected into 4 x 10° HEK293T cells using Lipofectamine 3000 (Thermo Fisher 
Scientific). After 4 h, the medium was replaced with fresh medium. Supernatants 
were collected 48 h after transfection and clarified by centrifugation at 3,000g, then 
passed through a 0.45-|um filter (Millipore). The filtered supernatants were stored 
at —80 °C in aliquots until use. To evaluate the incorporation of S proteins into the 
core of HIV virions, pseudoviruses in supernatant (20 ml) were concentrated by 
ultracentrifugation through a 2096 sucrose cushion (5 ml) at 80,000g for 90 min 
using a SW41 rotor (Beckman). Pelleted pseudoviruses were dissolved in 50 pl 
phosphate-buffered saline (PBS) and examined by electron microscopy. 
Pseudovirus infection. HeLa cells transiently expressing APN, ACE2 or DPP4 
were prepared using Lipofectamine 2000 (Thermo Fisher Scientific). Pseudoviruses 
prepared above were added to HeLa cells overexpressing APN, ACE2 or DPP4 
24 h after transfection. The unabsorbed viruses were removed and replaced with 
fresh medium at 3 h after infection. The infection was monitored by measuring 
the luciferase activity conferred by the reporter gene carried by the pseudovirus, 
using the Luciferase Assay System (Promega) as follows: cells were lysed 48 h after 
infection, and 20 ul of the lysates was taken for determining luciferase activity after 
the addition of 50 ul of luciferase substrate. 

Examination of known CoV receptors for SADS-CoV entry/infection. 
HeLa cells transiently expressing APN, ACE2 or DPP4 were prepared using 
Lipofectamine 2000 (Thermo Fisher Scientific) in a 96-well plate, with mock-trans- 
fected cells as controls. SADS-CoV grown in Vero cells was used to infect HeLa 
cells transiently expressing APN, ACE2 or DPP4. The inoculum was removed after 
1 h of absorption and washed twice with PBS and supplemented with medium. 
SARS-related-CoV WIV16’ and MERS-CoV HIV-pseudovirus were used as pos- 
itive control for human/swine ACE2 or human/swine DPP4, respectively. After 
24 h of infection, cells were washed with PBS and fixed with 496 formaldehyde in 
PBS (pH 7.4) for 20 min at room temperature. SARS-related-CoV WIV 16 rep- 
lication was detected using rabbit antibody against the SARS-related-CoV Rp3 
N protein (made in house, 1:100) followed by Cy3-conjugated goat anti-rabbit 
IgG (1:50, Proteintech)’. SADS-CoV replication was monitored using rabbit anti- 
body against the SADSr-CoV 3755 N protein (made in house, 1:50) followed by 
FITC-conjugated goat anti-rabbit IgG (1:50, Proteintech). Nuclei were stained with 
DAPI (Beyotime). Staining patterns were examined using confocal microscopy on 
a FV1200 microscope (Olympus). Infection of MERS-CoV HIV-pseudovirus was 
monitored by luciferase 48 h after infection. 

High-throughput sequencing, pathogen screening and genome assembly. Tissue 
from the small intestine of deceased pigs was homogenized and filtered through 
0.45-um filters before nucleic acid extraction and ribosomal RNA was depleted 
using the NEBNext rRNA Depletion Kit (New England Biolabs). Metagenomics 
analysis of both RNA and DNA viruses was performed. For RNA virus screening, 
the sequencing library was constructed using Ion Total RNA-Seq Kit v2 (Thermo 
Fisher Scientific). For DNA virus screening, NEBNext Fast DNA Fragmentation 
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& Library Prep Set for Ion Torrent (New England Biolabs) was used for library 
preparation. Both libraries were sequenced on an Ion S5 sequencer (Thermo 
Fisher Scientific). An analysis pipeline was applied to the sequencing data, which 
included the following analysis steps: (1) raw data quality filtering; (2) host genomic 
sequence filtering; (3) BLASTn search against the virus nucleotide database using 
BLAST; (4) BLASTx search against the virus protein database using DIAMOND 
v.0.9.0; (5) contig assembling and BLASTx search against the virus protein data- 
base. For whole viral genome sequencing, amplicon primers (provided upon 
request) were designed using the Thermo Fisher Scientific online tool with the 
HKU2-CoV and the SADS-CoV farm A genomes as references, and the sequenc- 
ing libraries were constructed using NEBNext Ultra II DNA Library Prep Kit for 
Illumina and sequenced on an MiSeq sequencer. PCR and Sanger sequencing was 
performed to fill gaps in the genome. Genome sequences were assembled using 
CLC Genomic Workbench v.9.0. 5'-RACE was performed to determine the 5’-end 
ofthe genomes using SMARTer RACE 5’/3’ Kit (Takara). Genomes were annotated 
using Clone Manager Professional Suite 8 (Sci-Ed Software). 








Phylogenetic analysis. SADS-CoV genome sequences and other representative 
coronavirus sequences (obtained from GenBank) were aligned using MAFFT 
v.7.221. Phylogenetic analyses with full-length genome, S gene and RdRp were 
performed using MrBayes v.3.2. Markov chain Monte Carlo was run for 20-50 
million steps using the GTR--G--I model (general time reversible model of nucleo- 
tide substitution with a proportion of invariant sites and y-distributed rates among 
sites). The first 1096 was removed as burn-in. The association between phylogenies 
and phenotypes (for example, host species and farms) was assessed by BaTS beta- 
build2, with the trees obtained in the previous step used as input. For SADS-CoVs, 
a median-joining network analysis was performed using PopART v.1.7, with e = 
0. Phylogenetic analysis of the 33 full-length SADS-CoV genome sequences was 
performed using RAxML v.8.2.11, with GTRGAMMA as the nucleotide substitu- 
tion model and 1,000 bootstrap replicates. The maximum likelihood tree was used 
to test the molecular clock using TempEst v.1.5. Potential genetic recombination 
events in our datasets were detected using RDP v.4.72. 

Animal infection studies. Experiments were carried out strictly in accordance 
with the recommendations of the Guide for the Care and Use of Laboratory 
Animals of the National Institutes of Health. The use of animals in this study 
was approved by the South China Agricultural University Committee of Animal 
Experiments (approval number 201004152). 

Two different animal challenge experiments were conducted. Pigs were used 
without gender preference. In the first experiment, which was conducted before 
the virus was isolated, we used three-day old specific pathogen-free (SPF) piglets 
of the same breeding line, cared for at a SPF facility, fed with colostrum (except 
one). These piglets were bred and reared to be free of PEDV, CSFV, SIV, PCV2 
and PPV infections, and were routinely tested for viral infections using PCR. We 
also conducted NGS to further confirm that these were animals were free of infec- 
tion of the above viruses before the animal experiment, and to demonstrate that 
the animals were free of SADS-CoV infection. The intestinal tissue samples from 
healthy and diseased animals (intestinal samples excised from euthanized piglets, 
then ground to make slurry for the inoculum and NGS was performed to confirm 
no other pig pathogens were found in the samples), were used to feed two groups 
of 5 (control) and 7 (infection) animals, respectively. For the second experiment, 
isolated SADS-CoV was used to infect healthy piglets from a farm in Guangdong, 
which had been free of diarrheal disease for a number of weeks. These piglets were 
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from the same breed as those on SADS-affected farms, to eliminate potential host 
factor differences and to more accurately reproduce the conditions that occurred 
during the outbreak in the region. Both groups of piglets were cared for at a known 
pig disease-free facility. Again, qPCR and NGS were used to make sure that there 
was no other known swine diarrhoea virus present in the virus inoculum or any 
of the experimental animals. Two groups (6 for each group) of three-day old pig- 
lets were inoculated with SADS-CoV culture supernatant or normal cell culture 
medium as control. NGS and qPCR were used to confirm that there were no other 
known swine pathogens in the inoculum. 

For both experiments, animals were recorded daily for signs of diseases, such 

as diarrhoea, weight loss and death. Faecal swabs were collected daily from all 
animals and screened for known swine diarrhoea viruses by qPCR. Weight loss was 
calculated as the percentage weight loss compared the original weight at day 0 with 
a threshold of 7596. It is important to point out that piglets when they are three 
days old tend to suffer from diarrhoea and weight loss when they are taken away 
from sows and the natural breast-feeding environment even without infection. At 
experimental endpoints, piglets were humanely euthanized and necropsies per- 
formed. Pictures were taken to record gross pathological changes to the intestines. 
Ileal, jejunal and duodenal tissues were taken from selected animals and stored at 
-80 °C for further analysis. 
Haematoxylin and eosin and immunohistochemistry analysis. Frozen (-80 °C) 
small intestinal tissues including duodenum, jejunum and ileum taken from the 
experimentally infected pigs were pre-frozen at -20 °C for 10 min. Tissues were 
then embedded in optimal cutting temperature (OCT) compound and cut into 
8-um sections using the Cryotome FSE machine (Thermo Fisher Scientific). 
Mounted microscope slides were fixed with paraformaldehyde and stained with 
haematoxylin and eosin for histopathological examination. 

For immunohistochemistry analysis, a rabbit antibody raised against the SADSr- 

CoV 3755 N protein was used for specific staining of SADS-CoV antigen. Slides 
were blocked by incubating with 10% goat serum (Beyotime) at 37 °C for 30 min, 
followed by overnight incubation at 4 °C with the rabbit anti-3755 N protein serum 
(1:1,000) and mouse anti-cytokeratin 84-18--19 monoclonal antibody (Abcam), 
diluted 1:100 in PBST buffer containing 596 goat serum. After washing, slides were 
then incubated for 50 min at room temperature with Cy3-conjugated goat-anti-rab- 
bit IgG (Proteintech) and FITC-conjugated goat-anti-mouse IgG (Proteintech), 
diluted 1:100 in PBST buffer containing 596 goat serum. Slides were stained with 
DAPI (Beyotime) and observed under a fluorescence microscope (Nikon). 
Reporting Summary. Further information on experimental design is available in 
the Nature Research Reporting Summary linked to this paper. 
Data availability. Sequence data that support the findings of this study have been 
deposited in GenBank with accession codes MF094681-MF094688, MF769416- 
MF769444, MF094697-MF094701, MF769406-MF769415 and MG557844. Raw 
sequencing data that support the findings of this study have been deposited in the 
Sequence Read Achieve (SRA) with accession codes SRR5991648, SRR5991649, 
SRR5991650, SRR5991651, SRR5991652, SRR5991654, SRR5991655, SRR5991656, 
SRR5991657, SRR5991658 and SRR5995595. 


2/. Harlow, E. & Lane, D. Antibodies: A Laboratory Manual (Cold Spring Harbor 
Laboratory Press, New York, 1988). 

28. Ren, W. et al. Difference in receptor usage between severe acute respiratory 
syndrome (SARS) coronavirus and SARS-like coronavirus of bat origin. J. Virol. 
82, 1899-1907 (2008). 
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Extended Data Fig. 1 | Map of outbreak locations and sampling sites SADS-CoV (sample 162140) originated in Conghua. The red flag marks 

in Guangdong province, China and the co-circulation of PEDV and Foshan city, the site of the SARS index case. b, Pooled intestinal samples (n 
SADS-CoV during the initial outbreak on farm A. a, SADS-affected — 5 or more biological independent samples) were collected at dates given 
farms are labelled (farms A-D) with blue swine silhouettes following the on the x axis from deceased piglets and analysed by qPCR. The viral load 
temporal sequence of the outbreaks. Bat sampling sites are indicated with for each piglet is shown as copy number per milligram of intestine tissue 
black bat silhouettes. The bat SADSr-CoV that is most closely related to (y axis). 


© 2018 Macmillan Publishers Limited, part of Springer Nature. All rights reserved. 


RESEARCH 








LETTER 








MF094681 SADS-CoV Farm A fi 
1 MF094682 SADS-CoV Farm B fi 
MF094683 SADS-CoV Farm C fi 
MF094684 SADS-CoV Farm D @& 
MF094687 SADSr-CoV 141388 «Wie 
MF094688 SADSr-CoV 162140 «i 
MF094686 SADSr-CoV 8495 «uie 
MF094685 SADSr-CoV 8462 «Wi 
1 NC 009988.1 Bat coronavirus HKU2 *ije 
NC 028811.1 BtMr-AlphaCoV SAX201 1 «we 
ır NC 002645.1 Human coronavirus 229E f 
:[- NC 028752.1 Camel alphacoronavirus isolate camel/Riyadh/Ry141/2015 Vfh 
i KY073747.1 229E-related bat coronavirus strain BEKY229E-1 «mfg 
' 7— NC 005831.2 Human Coronavirus NL63 k 
NC 032107.1 NL63-related bat coronavirus strain BEKYNL63-9a «ipt 
i NC 018871.1 Rousettus bat coronavirus HKU10 ij 
NC 028814.1 BtRf-AlphaCoV HuB2013 «wee 
Er NC 010437.1 Bat coronavirus 1A «ife 
1 F NC 010438.1 Bat coronavirus HKU8 «wee 
NC 028833.1 BtNv-AlphaCoV SC2013 «uie 
AF353511.1 PEDV strain CV777 fili 
KJ158152.1 PEDV strain AH-M ff 
JN547228.1 PEDV strain CH S fik 
1 JX647847.1 PEDV strain GD-1 fllik 
1 KJ960178.1 PEDV clone VN/JFP1013 1/2013/Vinh An/Vietnam fli 
‘I! KF761675.1 PEDV strain CH/YNKM-8/2013 fl 
KF272920.1 PEDV strain USA/Colorado/2013 fili 
KJ399978.1 PEDV strain OH851 fili 
KR061458.1 PEDV strain Italy/7239/2009 fli 
4 KU836638.1 PEDV isolate Yorkshire-2000 filik 
NC 009657.1 Scotophilus bat coronavirus 512 «wife 
NC 022103.1 Bat coronavirus CDPHE15/USA/2006 «wife 
1 ıı AJ271965.2 Transmissible gastroenteritis virus fili 
1 NC 028806.1 Swine enteric coronavirus strain Italy/213306/2009 fIlik 
4 E JQ404409.1 Canine coronavirus strain 1-71 Je 
NC 002306.3 Feline infectious peritonitis virus @ 
4 HM245926.1 Mink coronavirus strain WD1 133 -€ 
1 NC 023760.1 Mink coronavirus strain WD1127 -€ 
NC 030292.1 Ferret coronavirus isolate FRCoV-NL-2010 -€ 
NC 032730.1 Lucheng Rn rat coronavirus isolate Lucheng-19 4m. 
NC 034972.1 Coronavirus AcCoV-JC34 4m., 
NC 001846.1 Mouse hepatitis virus strain MHV-A59 C12 mutant 4a 
1 1. p NC 003045.1 Bovine coronavirus - 
1t KT368891.1 Camel coronavirus HKU23 isolate camel/Riyadh/Ry123/2015 he) 
DQ011855.1 Porcine hemagglutinating encephalomyelitis virus strain VW572 fik 
r NC 009019.1 Bat coronavirus HKU4-1 «x> 
1 i NC 009020.1 Bat coronavirus HKU5-1 «wipe 
7 KC869678.4 Coronavirus Neoromicia PML-PHE1/RSA/2011 «wee 
NC 019843.3 Middle East respiratory syndrome coronavirus $ 
NC 004718.3 SARS coronavirus $ 
C881005.1 Bat SARS-like coronavirus RSSHC014 «wipe 
F367457.1 Bat SARS-like coronavirus WIV1 «mig 
DQ022305.2 Bat SARS coronavirus HKU3-1 «wee 
NC 009021.1 Bat coronavirus HKU9-1 «wee 


E 


SADS cluster 


Alphacoronavirus 


Betacoronavirus 





0.3 


MF094681 SADS-CoV Farm A fiM 
MF094682 SADS-CoV Farm B fil 
MF094683 SADS-CoV Farm C fik 
MF094684 SADS-CoV Farm D 4l 
os] MF094687 SADSr-CoV 141388 pm 
MF094688 SADSr-CoV 162140 ie 
1 ,MF094685 SADSr-CoV 8462 «ue 
NC 009988.1 Bat coronavirus HKU2 «We 
MF094686 SADSr-CoV 8495 «wee 
,NC 002645.1 Human coronavirus 229E k 
[t NC 028752.1 Camel alphacoronavirus isolate camel/Riyadh/Ry141/2015 Mf 
. KY073747.1 229E-related bat coronavirus strain BEKY229E-1 «wipe 
NC 005831.2 Human Coronavirus NL63 f 
NC 032107.1 NL63-related bat coronavirus strain BEKYNL63-9a «wife 
NC 018871.1 Rousettus bat coronavirus HKU10 «wee 
, NC 028814.1 BtRf-AlphaCoV HuB2013 «wipe 
AF353511.1 PEDV strain CV777 fik 
KJ158152.1 PEDV strain AH-M iik 
1 JN547228.1 PEDV strain CHS fli 
JX647847.1 PEDV strain GD-1 
1 KJ960178.1 PEDV clone VN/JFP1013 1/2013/Vinh An/Vietnam fili 
KF272920.1 PEDV strain USA/Colorado/2013 fili 
KJ399978.1 PEDV strain OH851 #@& 
KF761675.1 PEDV strain CH/YNKM-8/2013 fik 
KR061458.1 PEDV strain Italy/7239/2009 4i 
4 KU836638.1 PEDV isolate Yorkshire-2000 fü 
NC 009657.1 Scotophilus bat coronavirus 512 aie 
NC 022103.1 Bat coronavirus CDPHE15/USA/2006  *W:ire 
- NC 010437.1 Bat coronavirus 1A *ir 
i NC 010438.1 Bat coronavirus HKU8 S% 
NC 028833.1 BtNv-AlphaCoV SC2013 = 
NC 028811.1 BtMr-AlphaCoV SAX2011 «wipe 
NC 032730.1 Lucheng Rn rat coronavirus isolate Lucheng-19 4m, 
NC 034972.1 Coronavirus AcCoV-JC34 4m, 
1 AJ271965.2 Transmissible gastroenteritis virus fik 
!|- NC 028806.1 Swine enteric coronavirus strain Italy/213306/2009 fI 
i 1| - JQ404409.1 Canine coronavirus strain 1-71 W€ 
NC 002306.3 Feline infectious peritonitis virus & 
,HM245926.1 Mink coronavirus strain WD1133 af 
4 - NC 023760.1 Mink coronavirus strain WD1127 .€ 
NC 030292.1 Ferret coronavirus isolate FRCoV-NL-2010 -€ 
; NC 001846.1 Mouse hepatitis virus strain MHV-A59 C12 mutant 4m. 
1 4 NC 003045.1 Bovine coronavirus FH 
1} KT368891.1 Camel coronavirus HKU23 isolate camel/Riyadh/Ry123/2015 wv 
DQ011855.1 Porcine hemagglutinating encephalomyelitis virus strain VW572 4@& 
: NC 009019.1 Bat coronavirus HKU4-1 «wipe 
1 ; NC 009020.1 Bat coronavirus HKU5-1 ije 
; KC869678.4 Coronavirus Neoromicia PML-PHE1/RSA/2011 «wie 
NC 019843.3 Middle East respiratory syndrome coronavirus $ 
if NC 004718.3 SARS coronavirus f 
KC881005.1 Bat SARS-like coronavirus RSSHCO14 «wir 
KF367457.1 Bat SARS-like coronavirus WIV1 «mpm 
DQ022305.2 Bat SARS coronavirus HKU3-1 «Wee 
NC 009021.1 Bat coronavirus HKU9-1 «wipe 


SADS cluster 


Alphacoronavirus 


Betacoronavirus 





0.3 


Extended Data Fig. 2 | Bayesian phylogenetic tree of the full-length 
genome and the ORF1a and ORFIb sequences of SADS-CoV and related 
coronaviruses. a, Bayesian phylogenetic tree of the full-length genome. 

b, Bayesian phylogenetic tree of the ORF1a and ORFIb sequences. Trees 
were constructed using MrBayes with the average standard deviation of 


split frequencies under 0.01. The host of each sequence is represented as 
a silhouette. Newly sequenced SADS-CoVs are highlighted in red, bat 
SADSr-CoVs are shown in blue and previously published sequences are 
shown in black. Scale bars, nucleotide substitutions per site. 
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Extended Data Fig. 3 | Phylogeny and haplotype network analyses of 
the 33 SADS-CoV strains from the four farms. a, Phylogenetic tree 
constructed using MrBayes. The GTR+#GAMMA model was applied 

and 20 million steps were run, with the first 10% removed as burn in. 
Viruses from different farms are labelled with different colours. Scale bar, 
nucleotide substitutions per site. b, Median-joining haplotype network 
constructed using ProART. In this analysis, € = 0 was used. The size of the 
circles represents the number of samples. The larger the circle, the more 


samples it includes. 
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Extended Data Fig. 5 | Isolation and antigenic characterization of rabbit serum raised against the recombinant SADSr-CoV N protein (red) 
SADS-CoV. a, b, Vero cells are shown 20 h after infection with mock (a) and DAPI (blue). The experiment was conducted independently three 


or SADS-CoV (b). c, d; Mock or SADS-CoV-infected samples stained with times with similar results. Scale bars, 100 jum. 
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Extended Data Table 1 | List of all known swine viruses tested by PCR at the beginning of the of SADS outbreak investigation on the four 
farms 








PEDV PDCoV TGEV RV PBV PSV SVA SIV NADC30 PRV FMDV CSFV PCV2 PCV3 APPV PPV Norovirus 


Farm A - - = : : = - E E M : i H 2 i ND 5 
Farm B - - ^ - - s j S i - E : j j ND H 
Farm C - - = . : = s F : - g E i 7 - 7 ND 
Farm D - - a = = : : = s : : 3 ; B E E ND 


Faeces, intestine or faecal swabs collected from January to April 2017 were tested. Sampling type and number of samples per farm were as follows. Farm A: 1 fecal sample, 20 intestinal sample and 6 
faecal swabs; farm B: 1 faecal sample and 15 intestinal samples; farm C: 2 intestinal sample and 1 faecal swab; farm D: 5 faecal sample and 1 faecal swab. The dash indicates a negative PCR result. 
ND, not determined. APPV, atypical porcine pestivirus; CSFV, classical swine fever virus; FMDV, foot and mouth disease virus; NADC30, porcine reproductive and respiratory syndrome virus, strain 
NADC30; PBV, porcine picobirnavirus; PCV2, porcine circovirus 2; PCV3, porcine circovirus 3; PDCoV, porcine deltacoronavirus; PEDV, porcine epidemic diarrhoea virus; PPV, porcine parvovirus; PRV, 
porcine pseudorabies virus; PSV, porcine sapelovirus; RV, porcine rotavirus; SIV, swine influenza virus; SVA, porcine senecavirus A; TGEV, porcine transmissible gastroenteritis virus. 
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Extended Data Table 2 | List of SADSr-CoVs detected in bats in Guangdong, China 


Sampling 


Time (Month-Year) 


Jun 13 


Jul 13 


Jul 13; May 14; 
Jun 15; Aug 16 


Sep 14; Jun 15; 
Aug 16 


Jun 13; Nov 13; 
Aug 14; Jun 15 


Jun 15 
Sep 14 


Location 


Yingde 


Yangshan 


Ruyuan 


Conghua 


Huidong 


Baoan 


Xiangzhou 


Bat Species 
Rhinolophus sinicus 
Pipistrellus abramus 
Myotis ricketti 
Pipistrellus abramus 
Hipposideros pratti 
Rhinolophus sinicus 
Rhinolophus affinis 
Rhinolophus macrotis 
Rhinolophus pusillus 
Rhinolophus rex 
Hipposideros pratti 
Rhinolophus sinicus 
Rhinolophus affinis 
Rhinolophus pusillus 
Hipposideros pomona 
Myotis ricketti 
Rhinolophus sinicus 
Rhinolophus affinis 
Rhinolophus macrotis 
Rhinolophus pusillus 
Hipposideros pomona 
Myotis ricketti 
Rhinolophus sinicus 
Rhinolophus pusillus 


Hipposideros pomona 
Total 


See Extended Data Fig. 1 for sampling sites in relation to SARS and SADS outbreak locations. 


PCR analysis 





LETTER Beye ea 








Fecal swabs sampled PCR Positive 


1 
8 
2 
1 
36 
2/ 
11 
3 
41 
9 
T 
70 
34 
11 
10 
1 
3/ 
59 
15 
1 
2 
84 
55 
28 
38 
591 


© 2018 Macmillan Publishers Limited, part of Springer Nature. All rights reserved. 


0 


BRIO O NMN NION WOO gio OJo oO 


ool|oeooeooNM 


cO 


58 (9.8%) 





asus uml LETTER 


Extended Data Table 3 | Test of SADS-CoV entry and infection in Hela cells expressing known coronavirus receptors 
HuAPN* HuACE2*  HuDPP4*  SwAPN* SwACE2*  SwDPP4* 








SADS-CoV : : : - - - 

SARS-related-CoV NA + NA NA + NA 
MERS-CoV'! NA NA + NA NA NA 
Expression! * (S-ta * (HIS-ta * (S-ta * (S-ta * (S-ta * (S-ta 





*Gene accession numbers for the genes used in this study: human APN, M22324.1; human ACE2, NM. 021804; human DPP4, NM. 001935.3; SwAPN (swine APN), NM 214277.1; SwACE2 (swine 
ACE2), NM. 001116542.1; SwDPP4 (swine DPP4), NM 214257.1. 

TFor MERS-CoV infection, HIV-pseudovirus was used. 

+Expression of APN, ACE2 and DPP4 was confirmed by antibodies against fused tags. 
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Extended Data Table 4 | Experimental infection of SPF piglets using intestine tissue homogenate 








a 
Group Animal Age Inoculum SADS-CoV Inoculum Inoculation Data recorded on day one and (day two) post challenge 
Number (days) material titer volume route Death Weight Watery SADS-CoV PEDV/PDCoV/RV 
(copy/ml) loss diarrhea +ve +ve 
Infected 7 3 PCR positive intestine slurry 1.55x10*6 4 ml Oral + milk 0/7 (3/7) 4/7 (5/7) 5/7 (7/7) 6/7 (7/7) 0/7 (0/7) 
Control 5 3 PCR negative intestine slurry 0 4 ml Oral + milk 0/5 (1/5) 1/5 (3/5) 0/5 (1/5) 0/5 (0/5) 0/5 (0/5) 
b 
Days post p : : " - M : " - t . aet : ai 
Group challenge Piglet-I1 Piglet-I2 Piglet-I3 Piglet-l4 Piglet-I5 Piglet-l6 Piglet-I7 
0 0.565 0.66 0.6 0.68 0.49 0.57 0.62 
Infected 1 0.555 0.635 0.685 0.715 0.4 0.475 0.565 
2 0.51 0.52 0.665 0.785 


Piglet-C1* ^ Piglet-C2*  Piglet-C3*  Piglet-C4" Piglet-C5* 


0 0.67 0.59 0.5 0.53 0.525 
Control 1 0.765 0.53 0.49 0.51 0.535 
2 0.765 0.53 0.575 0.505 


Experimental details can be found in the Methods. a, Animals were recorded every day for signs of disease, including weight loss, diarrhoea and death. PCR on DNA from faecal swabs was carried out 
to monitor the presence of SADS-CoV or other pig viruses. b, Daily body weight record of all piglets. Weights are in kg. 
*Euthanized on the indicated day for further analysis. 


TAnimal died during the experiment. 
tiThe only animal that did not receive colostrum in this experiment due to shortage in supply. 
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Extended Data Table 5 | Experimental animal infection of farm piglets using cultured SADS-CoV 








a 
~ Group Anima = Age  lnocuum  SADS-CoV Inoculum Inoculation — Data recorded on day two and (day four) post challenge ^ ^ — 
Number (days) material titer volume route Death Weight Watery SADS-CoV PEDV/PDCoV/RV 
(TCID;,/ml) loss diarrhea +ve +ve 
Infected 6 3 Cultured SADS-CoV 10*6.625 6 ml Oral * milk 1/6 (3/6) 4/6 (6/6) 6/6 (6/6) 6/6 (6/6) 0/6 (0/6) 
Control 6 3 Mock culture supernatant 0 6 ml Oral * milk 0/6 (0/6) 3/6 (3/6) 5/6 (3/6) 0/6 (0/6) 0/6 (0/6) 
b 
Group Days post challenge Piglet-I1' Piglet-I2' Piglet-13* Piglet-I4* Piglet-15* Piglet-I6' 
0 1.5 1.54 2.32 1.92 1.54 2.165 
1 1.41 1.575 2.58 1.885 1.46 2.08 
Infected 2 1.23 1.39 2.615 1.73 1.54 1.365 
3 2.115 1.54 1.335 1.725 
4 1.505 
Piglet-C1* ^ Piglet-C2* Piglet-C3* ^ Piglet-C4* ^ Piglet-C5*  Piglet-C6* 
0 1.955 2.055 2.8 1.835 1.835 1.83 
1 1.765 1.955 1.9 1.68 1.645 1.93 
Control 2 2.12 1.675 1.93 1.515 1.9 
3 2.25 1.69 2.18 1.66 2.38 
4 2.27 1.555 2.58 


Experimental details can be found in the Methods. a, Animals were recorded every day for signs of disease, including weight loss, diarrhoea and death. PCR on DNA from faecal swabs was carried out 
to monitor the presence of SADS-CoV or other pig viruses. b, Daily body weight record of all piglets. Weights are in kg. 

*Euthanized on the indicated day for further analysis. 

TAnimal died during the experiment. 
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Nature Research wishes to improve the reproducibility of the work that we publish. This form is intended for publication with all accepted life 
science papers and provides structure for consistency and transparency in reporting. Every life science submission will use this form; some list 
items might not apply to an individual manuscript, but all fields must be completed for clarity. 


For further information on the points included in this form, see Reporting Life Sciences Research. For further information on Nature Research 
policies, including our data availability policy, see Authors & Referees and the Editorial Policy Checklist. 


>» Experimental design 
1. Sample size 


Describe how sample size was determined. For the three main figures: figure 1 calculated all sick pigs or used more than five 
samples in epidemiology, or used three animals in tissue distribution (meets the 
minimal statistical requirements). Figure 2 used most of the representative CoV 
genomes thus should be adequate. For all other tables or figures that sample size 
involved, we used more than three samples per group. For animal experiments, we 
used at least five animal per group. 


2. Data exclusions 


Describe any data exclusions. No data exclusion. 

3. Replication 
Describe whether the experimental findings were As epidemiology study, we presented all results including positive or negative here. 
reliably reproduced. The authors guarantee the findings are reliably reproducible. At least three 


independent experiments were performed, which was stated in the text. 


4. Randomization 


Describe how samples/organisms/participants were Animals were randomly assigned to groups prior to any experimentation. 
allocated into experimental groups. 

5. Blinding 
Describe whether the investigators were blinded to SADS-CoV histology was performed in a blinded manner. 


group allocation during data collection and/or analysis. 


Note: all studies involving animals and/or human research participants must disclose whether blinding and randomization were used. 


6. Statistical parameters 


For all figures and tables that use statistical methods, confirm that the following items are present in relevant figure legends (or in the 
Methods section if additional space is needed). 


n/a | Confirmed 


LI x] The exact sample size (n) for each experimental group/condition, given as a discrete number and unit of measurement (animals, litters, cultures, etc.) 


E x< A description of how samples were collected, noting whether measurements were taken from distinct samples or whether the same 
sample was measured repeatedly 


||| ><] A statement indicating how many times each experiment was replicated 


d E The statistical test(s) used and whether they are one- or two-sided (note: only common tests should be described solely by name; more 
complex techniques should be described in the Methods section) 


x] L| A description of any assumptions or corrections, such as an adjustment for multiple comparisons 
Dx] [| The test results (e.g. P values) given as exact values whenever possible and with confidence intervals noted 
| | A clear description of statistics including central tendency (e.g. median, mean) and variation (e.g. standard deviation, interquartile range) 


| | Clearly defined error bars 


See the web collection on statistics for biologists for further resources and guidance. 
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> Software 


Policy information about availability of computer code 


7. Software 


Describe the software used to analyze the data in this BLAST+ v2.2.3, CLC Genomic Workbench v9.0, Clone Manager v8, MAFFT v7.221, 
study. MrBayes v3.2, DIAMOND v0.9.0, BaTS beta-build2, PopART v1.7, RAxML v8.2.11, 
TempEst v1.5, RDP v4.72. 


For manuscripts utilizing custom algorithms or software that are central to the paper but not yet described in the published literature, software must be made 
available to editors and reviewers upon request. We strongly encourage code deposition in a community repository (e.g. GitHub). Nature Methods guidance for 
providing algorithms and software for publication provides further information on this topic. 


» Materials and reagents 


Policy information about availability of materials 


8. Materials availability 


Indicate whether there are restrictions on availability of | There is no restriction to material availability. 
unique materials or if these materials are only available 
for distribution by a for-profit company. 
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Describe the antibodies used and how they were validated 1, rabbit anti-HKU2-NP polyclonal antibody, made by ourselves, validated by 

for use in the system under study (i.e. assay and species). immunogen in a WB (titer 1: 10000); 2, anti-HIS tag monoclonal antibody 
(Proteintect Group), validated in a WB (titer 1: 1000); 3, anti-S tag monoclonal 
antibody, made by ourselves, validated in a WB (titer 1: 10000); 4, cyanin 3-labeled 
goat anti-mouse/rabbit IgG (Proteintech Group), validated in IFA (titer 1: 1000); 5, 
mouse anti-FLAG tag antibody (Thermo Fisher Scientific), validated in a WB (titer 1: 
1000); 6, mouse anti-Cytokeratin 8418419 mAb (Abcam), validated in IHC (1:100); 
7, FITC conjugated goat-anti-rabbit IgG (Proteintech), validated in IHC (1:100); 


10. Eukaryotic cell lines 


a. State the source of each eukaryotic cell line used. 1, African green monkey origin, Vero from ATCC; 2, bat origin Rhinolophus sinicus 
(made by ourselves), Kidney primary RsKi9409, lung primary RsLu4323, lung 
immortalized RsLuT, brain immortalized RsBrT and heart immortalized RsHeT; all 
bats were made in house; 3, Swine cells: intestinal IPEC and SIEC, kidney PK15, LLC- 
PK1 and IBRS, testes cell ST; all swine cells were from ATCC; 4, human cells: Hela 
and HEK293T were from ATCC. 


b. Describe the method of cell line authentication used. | All monkey and human cells were from ATCC with authentication. Swine cells 
(commercially available) were gifts of collaborators and were originally from ATCC 
with authentication. They were authentication by microscope observation during 
culture. Bat cells made by ourselves were from organ or cultured and 
immortalized. We guarantee they were from the organs described but there was 
no further authentication. 


c. Report whether the cell lines were tested for We confirm that all cells were tested as mycoplasma negative. 
mycoplasma contamination. 


d. If any of the cell lines used are listed in the database None of the cell lines used are listed in the ICLAC database. 
of commonly misidentified cell lines maintained by 
ICLAC, provide a scientific rationale for their use. 


» Animals and human research participants 


Policy information about studies involving animals; when reporting animal research, follow the ARRIVE guidelines 


11. Description of research animals 


Provide details on animals and/or animal-derived Swine used in animal infection study aged between 2-4 days. The first experiment 

materials used in the study. used healthy Chinese Bamaxiang SPF piglets that were cultured free of SADS-CoV 
or other known swine disease agents. The second experiment used healthy duroc- 
landrace-yokshire piglets (not SPF) that were not affected by SADS-CoV before. No 
gender preference when choose the animal. Piglets were from same breed and at 
same age and were randomly assigned into groups for the experiments. 
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Policy information about studies involving human research participants 


12. Description of human research participants 


Describe the covariate-relevant population Pig farm workers were bleed for testing possible spillover of SADS-CoV. These 
characteristics of the human research participants. workers are also adult male who had close contact with sick pigs. Non of them had 
clinical signs of diseases during sampling. 
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